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Abstract. Cosmic ray (CR) transport and acceleration is essential for many as- 
trophysical problems, e.g., CMB foreground, ionization of molecular clouds and all 
high energy phenomena. Recent advances in MHD turbulence call for revisions in the 
paradigm of cosmic ray transport. We use the models of magnetohydrodynamic turbu- 
lence that were tested in numerical simulation, in which turbulence is injected at large 
scale and cascades to to small scales. We shall address the issue of the transport of 
CRs, both parallel and perpendicular to the magnetic field and show that the issue of 
cosmic ray subdiffusion is only important for restricted cases when the ambient tur- 
bulence is far from that suggested by numerical simulations. Moreover, on scales less 
than injection scale of turbulence, CRs's transport becomes super-diffusive. We also 
shall discuss the nonlinear growth of kinetic gyroresonance instability of cosmic rays 
induced by large scale compressible turbulence. This gyroresonance of cosmic rays 
on turbulence is demonstrated an important scattering mechanism in addition to direct 
interaction with the compressible turbulence. The feedback of the instability on large 
scale turbulence cannot be neglected, and should be included in future simulations. 



1. Introduction 

The propagation and acceleration of cosmic rays (CRs) is governed by their interactions 
with magnetic fields. Astrophysical magnetic fields are turbulent and, therefore, the 
resonant and non-resonant (e.g. transient time damping, or TTD) interaction of cosmic 
rays with MHD t urbulence is the acce pted principal mechanism to scatter and isotropize 
cosmic rays (see ISchlickeiseil l2002h . In addition, efficient scattering is essential for 



the acceleration of cosmic rays. For instance, scattering of cosmic ra ys back into th e 



shock is a vital component of the first order Fermi acceleration (see lLongairill997n . 
At the same time, stochastic acceleration by turbulence is entirely based on scattering. 
The dynamics of cosmic rays in MHD turbulence holds the key to all high energy 
astrophysics and related problems. 

We live in an exiting era when we are starting to test fundamental processes taking 
place at the beginning of the Universe, at the event horizon of black holes, when the 
nature of dark matter and dark energy is being probed etc. In the mean time, with the 
launching of the new facilities like Fermi, we have much more observational data avail- 
able than ever before. Using computers many researchers make sophisticated complex 
models to confront the observations in unprecedented details. This makes it urgent that 
we understand the key physical processes underlying astrophysical phenomena, can 
parameterize them and, if necessary, use as a subgrid input in our computer models. 
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At present, the propagation of the CRs is an advanced theory, which makes use 
both of analytical studies and numerical simulations. However, these advances have 
been done within the turbulence paradigm which is being changed by the current re- 
search in the field. Instead of the empirical 2D+slab model of turbulence, numerical 
simulations suggest anisotropic Alfvenic modes (an analog of 2D, but not an exact one, 
as the anisotropv chan ges with the scale involved) + fast modes or/and slab modes 
(ICho & Lazarianll2002h . 

At the same time, one should not disregard the possibilities of generation of addi- 
tional perturbations on small scales by CR themselves. For in stance, the slab Alfvenic 
pertu rbation can be created, e.g., via streaming instability (see IWentzelll 1 974l : ICesarsky 



198(]h . These perturbations are present for a range of CRs energies (e.g., < lOOGeV 



in interstellar medium for the streaming instabihty) owing t o non-linear damping aris- 
ing from ambient turbulen ce (lYan & LazarianI |2002|. 12004 , henceforth YL02, YL04, 
Farmer & GoldreichI |2004) . Instabilities induced by anisotropic distribution o f CRs 



were also suggested as a possibility to scatter CRs dLercha . 1 96?! : iMekosd 1 1 974) . Gy- 
roresonance instability induced by large scale turbulence compressions, ca n be an im- 



portant feedbac k proce sses that creates slab modes to efficiently scatter CRs (iLazarian & Beresnyak 
l2"006 : ,Yan & Lazarianll201lh . 

Propagation of CRs perpendicular to mean magnetic field is another important 
problem for which one needs to take into account both large and small scale interac- 
tions in tested models of turbulence. Indeed, if one takes only the diffusion along the 
magn etic field line and fi eld line random walk (FLRW Jokipii 1966; Jokipii & Parked 
ll969l : iForman et al.lll974l ). compound (or subdiffusion) would arise. Whethe r the subd- 



iffusion is realistic in fact depends on the models of turbulence chosen ( Yan & Lazarian 



I2OO8). In this paper we again seek the answer to this questions within domain of nu- 
merically tested models of MHD turbulence. 

In what follows, we discuss the cosmic ray transport in large scale turbulence in 
§2. In §3, we shall study the perpendicular transport of cosmic rays on both large and 
small scales. We shall also discuss the issue of super-diffusion and the applicability of 
sub-diffusion. In §4, we address the issue of gyroresonance instability of CRs and its 
feedback on large scale compressible turbulence. Summary is provided in §5. 



2. Cosmic ray scattering by compressible MHD turbulence 

As we mentioned earlier, numerical simulations of MHD turbulence supported the 
GS95 model of turbulence, which does not have the "slab" Alfvenic modes that pro- 
duced most of the scattering in the earlier models of CR propagation. Can the turbu- 
lence that does not appeal to CRs bac k-reaction (see §4) produce ef ficient sca ttering? 
In the models of ISM turbulence ([Armstrong et al.ll995. : .McKee & Ostriker„2007.) . 



where the injection happens at large scale, fast modes were identified as a scatter- 
ing agent for cosmic rays in interstellar medium 

(S 

& LazarianI 12002. 20041). These 



works made use of the quantitative description of turbulence obtained in C ho & LazarianI 
(I2OO2) to calculate the scattering rate of cosmic rays. For instance, the scattering rate 
of relativistic protons by Alfvenic turbulence was shown to be nearly 10^" times lower 
than the generally accepted estimates obtained assuming the Kolmogorov scaling of 
turbulence (see F ig. [J). Although this estimate is 10^ times larger than that obtained by 
IChandranI (120001) . who employed GS95 ideas of anisotropy, but lacked the quantitative 
description of the eddies, it is clear that for most interstellar circumstances the Alfvenic 
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scattering is suppressed. The low efficiency of scattering by Alfvenic modes arise from 
higli anisotropy of the modes at the scale s of cosm ic ray gyroradius. 

Fast modes, however, are isotropic (iCho & Lazarian 2002 .). Indeed they a r e sub- 
ject to both collisional and collisionless damping. The studies in lYan & LazarianI (12002 , 



20041) demonstrated, nevertheless, that the scattering by fast modes dominates in most 



cases in spite of the damping. More recent studies of cosmic ray propagation and ac- 




Figure 1. Left: rate of CR scattering by Alfven waves versus CR energy. The 
lines at the top of the figure are the accepted e stimates obtained for Kolmogorov 
turbulence. The dotted curve is from' Chandrad fcOOO ). The analytical calculations 
are given by the solid line with our numerical calculations given by crosses. Middle: 
dimensionless CR scattering coefficient D^i^/£l for the case of /j = 0.71 vs Larmor 
radius expressed in cube size units (solid line). We suppose that it is dominated 
by magnetic bottles formed by slow mode, this is why (dimensional scattering 
frequency) is almost constant, i.e. independent on particle's energy. For comparison, 
we plot various theoretical predictions: QLT pre diction for Alfven and slow mode 
(dashed) (from iBeresnyak. Yan. & Lazarian 2011); Right: mean free path of CRs 
vs. their energy. The scattering is dominated by direct interaction by fast modes for 
CRs with energies > a few tens GeV and through gyroresonance instability due to 
turbulence compression for lower energy CRs. 



celeration that explicitly appeal to th e effect of the fast r nodes i nclude 'Cassano & Brunettil 
(12005.) : .Brunetti & Lazigi^ (l2007l) : I Yan & LazarianI (l2008l) : lYan et al. (2Q0i). Inci- 
dentally, fast modes h ave been also identified as primary agent s for the acceleration of 
charged dust particles lYan & LazitfianI (l2003l) : lYan et all (|2004|) . 



3. Perpendicular transport 

In this section we deal with the diffusion perpendicular to mean magnetic field. Both 
observation of Galactic CRs and solar wind indicate that the diffusion of CRs perpen- 
dicular to magnetic fie ld is comparable to parallel diffusion (.Giacalone & Jokipii 19991 : 
Maclennan et al.ll200lh . Why is that? 

The perpendicular transport is slow if particles are restricted to the magnetic field 
lines and the transport is solely due to the random walk of field line wandering (see 



' On tlie basis of weak turbulence theory, IChandrai] ( l2005h has argued that high-frequency fast waves, 
which move mostly parallel to magnetic field, generate Alfven waves also moving mostly parallel to 
magnetic field. We expect that the scattering by thus generated Alfven modes to be similar to the scat- 
tering by the fast mod es created by them. Therefore we expect that the simplified approach adopted in 
lYan & LazarianI ( 120041) and the papers that followed to hold. 
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K6ta& Jokipiill2000h . In the three-dime nsional turbulence, field lines are diverging 



away due t o shearing by Alfven modes (see 'Lazarian & Vishniac"l999"; Narayan & Medvedev 



20011 : iLazar i an 2006). Since the Larmor radii of CRs are much larger than the minimum 



scale of eddies l±^min, field hnes within the CR Larmor orbit are effectively diverging 
away owing to shear by Alfvenic turbulence. The cross-field transport thus results from 
the deviations of field lines at small scales, as well as field line random walk at large 
scale (> min[L/M^,L]), where L,Ma are the injection scale and Alfvenic Mach number 
of the turbulence, respectively. 

3.1. Perpendicular dilTusion on large scale 

High Ma turbulence: High Ma turbulence corresponds to the field that is easily bended 
by hydrodynamic motions at the injection scale as the hydro energy at the injection 
scale is much larger than the magnetic energy, i.e. pV^ » B^. In this case magnetic 
field beco mes dynam ically important on a much smaller scale, i.e. the scale I a = L/M^ 



LazarianI (see 2006). If the parallel mean free path of CRs <C U, the stiffness of 
B field is negligible so that the perpendicular diffusion coefficient is the same as the 
parallel one, i.e., D± = D||. If ^ 1a, the diffusion is controlled by the straightness of 
the field lines, and D± = f» 1 /31av, Ma > 1, X^\> Ia - The diffusion is isotropic if 
scales larger than 1a are concerned. In the opposite limit X|| < 1a, naturally, a result for 
isotropic turbulence, namely, D± =D^\ 1/3X||V holds. 

Low Ma turbulence: For strong magnetic field, i.e. the field that cannot be easily 
bended at the turbulence injection scale, individual magnetic field lines are aligned 
with the mean magnetic field. The diffusion in this case is anisotropic. If turbulence 
is injected at scale L it stays weak for the scales larger than LM\ and it is strong at 
smaller scales. Consider first the case of X\\ > L. The time of the individual step is 
L/v||, then ~ l/3LvM^, Ma < 1, X|| > L. This is similar to the case discussed in 
the FLRW model (Jokipii 1966). However, we obtain the dependence of M^ instead of 
their scaling. In the opposite case of < L, the perpendicular diffusion coefficient 
is D± ^ ^ll^A' which coincide s with the resul t obtained for the diffusion of thermal 
electrons in magnetized plasma (iLazarianll2006h . 



3.2. Superdiffusion on small scales 

The diffusion of CR on the scales ^ L is different and it is determined by how fast 
field lines aie diverging away fro m each other. The mean deviation of a field in a dis- 
tance 5x is proportional to [5z]^/^ dLazarian & Vishniac 119991 : iLazarianl I2006h. same as 



Richardson diffusion in the case of hy drodynamic turbulence (see lEyink et all 1201 Ih . 



Following the argument, we showed in I Yan & LazarianI ( 2008*) that the cosmic ray per- 



pendicular transport is superdiffusive. The reason is that there is no random walk on 
small scales up to the injection scale of strong MHD turbulence {LM^ for Ma < 1 and 
1a for Ma > 1). 

3.3. Is there subdiffusion? 

The diffusion coefficient D^^M'^ we obtained in the case of Ma < 1, means that the trans- 
port perpendicular to the dynamically strong magnetic field is a diffusion, rather than 
subdiffusion, as it was st ated in a number of re cent papers. This is also supported by 
test particle simulations dfieresnyak et al.ll201 ij) . Let us clarify this point by obtaining 



the necessary conditions for the subdiffusion to take place. 
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The major implicit assumption in subdiffusion (or compound diffusion) is that the 
particles trace back their trajectories in x direction on the scale 6z. When is it possible 
to talk about tracing particle trajectories back? In the case of random motions at a 
single scale only, the distance over which t he particle trajectories get uncorrected is 
given by the Rech ester & RosenbluthI (Il978h model. Assuming that the damping scale 
of the turbulence is larger that the CR Larmor radius, this model, when gene ralized 
to anisotropic turbulence provides (iNarayan & Medvedevll200ll : lLazarianll2006l) Lrr = 
l\\,min^^{t±,min/fLar) whcrc /[| „„„ is the parallel scale of the cut-off of turbulent motions, 
l±,min is the corresponding perpendicular scale, r/^^ is the CR Larmor radius. The 
assumption of r^a^ < l±,min can be valid, for instance, for Alfvenic motions in partially 
ionized gas. However, it is easy to see that, even in this case, the corresponding scale is 
rather small and therefore subdiffusion is not applicable for the transport of particles in 
Alfvenic turbulence over scales » l\\.min- 

If riar > l±,min, as it is a usual case for Alfven motions in the phase of ISM with 
the ionization larger than « 93 %, where the A lfvenic motio ns go to the the rmal parti- 
cle gyroradius (see estimates in lLithwick & G oldreich 200l1: lLazarian et aL 112004), the 
subdiffusion of CR is not an applicable concept for Alfvenic turbulence. This does 
not preclude subdiffusion from taking pla ce in particul ar models of magnetic perturba- 
tions, e.g. in the slab model considered in lKota & Jokip ii (200 0|), but we believe in the 
omnipresence of Alfvenic turbulence in interstellar gas (see lArmstrong et aTll 19951) . 



4. Gyroresonance Instability of CRs in Compressible Turbulence 

Until recently, test particle approximation was assumed in most of earlier studies and no 
feedback of CRs is included apart from the streaming instability. Turbulence cascade 
is established from large scales and no feedback of CRs is included. This may not 
reflect the reality as we know the ener gy of CRs is comparabl e to that in turbulence 
and m agnetic field (see iKulsrudI l2005h . It was suggested by iLazarian & Beresnyak 



(12006") that the gyroresonance instability of CRs can drain energy from the large scale 
turbulence and cause instability on small scales by the turbulence compression induced 
anisotropy on CRs. And the wave generate d on the scales, in turn, provides additional 
scattering to CRs. In lYan & LazarianI (I2O110 . we provided quantitative studies based on 



the nonlinear theory of the growth of the instability to take into account more accurately 
the feedback of the growing waves on the distributions of CRs. 

In the presence of background compressible turbulence, the CR distribution is 
bound to be anisotropic because of the conservation of the first adiabatic invariant 
fj = v\/B. Such anisotropic distribution is subjected to various instabilities. Waves 
are generated through the instabilities, enhancing the scattering rates of the particles, 
their distribution will be relaxed to the state of marginal state of instability even in the 
collisionless environment. While the hydrodynamic instability requires certain thresh- 
old, the kinetic instability can grow very fast with small deviations from isotropy. Here, 
we focus on the gyroresonance instability. Both the qualitative and quantitative stud- 
Yan & LazarianI (1201 ih show that the isotropization rate is roughly xJcatt ~ ^-^^ 



les m 



where Tcr,^N are the instability growth rate and the wave energy normalized by mag- 
netic energy, respectively. Pes is the ratio of CR pressure to magnetic pressure, A is the 
degree of anisotropy of the CR momentum distribution. 
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By balancing the rate of decrease in anisotropy due to scattering and the growth 
due to compression, one can get 



where va is the Alfven speed, co, 5v are the wave frequency and amplitude at the scale 
that effectively compresses the magnetic field and create anisotropy in CRs' distribution 
(lYan&Lazarian.2011,) . 

4.1. bottle-neck for the growth of the instability and feedback on turbulence 

The creation of the slab waves through the CR resonant instability is another channel 
to drain the energy of large scale turbulence. This process, on one hand, can damp the 
turbulence. On the other hand, it means that the growth rate is limited by the turbulence 
cascade. The energy growth rate cannot be larger than the turbulence energy cascading 
rate, which is 1 /2pV^ /va/L for fast modes in low p medium and Pv\/Ia for slow modes 
in high P medium. This places a constraint on the growth, thus the upper limit of wave 
energy is given by 



where y is the Lorentz factor and L,- ~ 6.4 x 10"'(B/5^G)(10" 10cm^/?icr)pc. The 
growth is induced by the compression at scales < Xcr- Therefore, in the case that 
TgrE reaches the energy cascading rate, fast modes are damped at the corresponding 
maximum turbulence pumping scale Xfb = rp/zp^. If X/^ is larger than the original 
damping scale Ic, then there is a feedback on the large scale compressible turbulence. 
This shows that test particle approach is not adequate and feedback should be included 
in future simulations. 

5. Summary 

• Compressible fast modes are most important for CR scattering. CR transport 
therefore varies from place to place. 

• Large scale mirror is essential for pitch angle scattering ( including 90 degree). 

• Subdiffusion does not happen in 3D turbulence . 

• Our results are tested using input from turbulence simulations. 

• Small scale slab waves are generated in compressible turbulence by gyrores- 
onance instability, dominating the scattering of low energy CRs (< lOOGeV). 
Feedback of CRs on turbulence should be included in future simulations. 



TgrVA 



(1) 




P<1 
P>1, 
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